From Gauging Nonrelativistic Translations to 

N-Body Dynamics 

J. Lukierski* 

O" 

Institute for Theoretical Physics, University of Wroclaw, 
(N ! pi. Maxa Borna 9, 50-204 Wroclaw, Poland 



> 
(N 
O 



Oh 
43 



X 



P.C. Stichelt 

An der Krebskuhle 21 
D-33619 Bielefeld 



■ W.J. Zakrzewski* 

O 

o 
o 



Department of Mathematical Sciences, University of Durham, 

Durham DH1 3LE, UK 



Abstract 

We consider the gauging of space translations with time-dependent gauge 
functions. Using fixed time gauge of relativistic theory, we consider the gauge- 
invariant model describing the motion of nonrelativistic particles. When we use 
gauge-invariant nonrelativistic velocities as independent variables the transla- 
tion gauge fields enter the equations through a d x (d + 1) matrix of vielbein 
fields and their Abelian field strengths, which can be identified with the tor- 
sion tensors of teleparallel formulation of relativity theory. We consider the 
planar case (d = 2) in some detail, with the assumption that the action for 
the dreibein fields is given by the translational Chern-Simons term. We fix 
the asymptotic transformations in such a way that the space part of the met- 
ric becomes asymptotically Euclidean. The residual symmetries are (local in 
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time) translations and rigid rotations. We describe the effective interaction 
of the d = 2 iV-particle problem and discuss its classical solution for N = 2. 
The phase space Hamiltonian H describing two-body interactions satisfies a 
nonlinear equation H = 7i(x,p; H) which implies, after quantization, a non- 
standard form of the Schrodinger equation with energy dependent fractional 
angular momentum eigenvalues. Quantum solutions of the two-body problem 
are discussed. The bound states with discrete energy levels correspond to a 
confined classical motion (for the planar distance between two particles r < ro) 
and the scattering states with continuum energy correspond to the classical 
motion for r > tq. We extend our considerations by introducing an external 
constant magnetic field and, for N = 2, provide the classical and quantum 
solutions in the confined and unconfined regimes. 

1 Introduction 

Our aim here is to discuss theories invariant under local time-dependent nonrelativis- 
tic translations (x — (x± . . .Xd))'- 

Xi = x'i(Z,t) (1.1) 

supplemented by global space rotations (x- = R/xj) and global time translations 

t' = t + a. (1.2) 

The usual approach to local coordinate invariance in nonrelativistic theory is to con- 
sider the limit c — > oo of a relativistic generally covariant theory. In particular it 
should be pointed out that the nonrelativistic limit of the Einstein action coupled to 
a relativistic point particle had been explicitly performed by Lusanna et al [1]. In 
this paper, following the earlier treatments by one of the present authors in the case 
of one-dimensional model [2,3], we shall impose nonrelativistic framework directly by 
constructing actions covariant under (1.1-2). 

We shall consider the problem here by using the vielbein formulation of (d + 
1) -dimensional relativistic gravity, with d vectors E~ (a = 1, . . . , d; // = 0, 1, . . . d) 
describing translational gauge fieldsQ obtained from a relativistic (d + 1) x (d + 1) 
vielbein Eji (p = 0, 1, . . . d) by fixing completely E^. We will find that in the equation 
of motion for gauge-invariant nonrelativistic velocities 

e = E%x t + Ei, (1.3) 

the derivative of the vielbein field will enter only through its Abelian field strength 

V H = - d v El (1.4) 

1 This idea goes back to the papers by Cho [4], Hayashi and Shirafuji [5,6]; for a review see [7]. 
Such framework leads to the so-called teleparallel formulation of relativity (see e.g. [5-10]), with 
vanishing curvature and nonvanishing torsion. 
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thus leading to the interpretation in terms of components of the torsion tensor. 

The aim of this paper is to consider the dynamical consequences of the coupling 
of nonrelativistic particles in two space dimensions to nonstandard D = 2 + 1 grav- 
ity action. It has been found in analogous one-dimensional model [2,3] that gauging 
of nonrelativistic translations with Maxwell-like field action quadratic in the torsion 
tensor leads to confinement via the geometric bag formation. Further one can pos- 
tulateP] that such a confinement can occur on a new level of microscopic theories of 
fundamental interactions and does not appear in macroscopic "physical" gravity. The 
main result of the present paper is to showQ that the geometric bag solutions which 
may describe confinement occur also in the planar case. 

Our paper is organized as follows. Firstly, in Sect. 2 we shall consider d- dimensional 
nonrelativistic point particles, by using the Lagrangian framework covariant under 
(1.1). After considering in more detail the reparametrization-invariant nonrelativis- 
tic particle dynamics for d — 1, 2 and 3, in Sect. 3 we shall consider the field actions 
for vielbein fields, using known results from the teleparallel formulation of relativistic 
gravity in 1 + 1, 2 + 1 and 3 + 1 dimensions. We shall find that the d = 1 action 
proposed in [2] is obtained from the well-known action in D = 1 + 1 Einstein-Cartan 
gravity [14]; for D = 2 + 1 we will consider the field action given by the so-called 
translational Chern-Simons term [15-17]. 

In Sect. 4 we solve the field equations having chosen the appropriate boundary 
conditions (gauge fixing). The general considerations of the iV-particle classical par- 
ticle dynamics in Sect. 5 is specialized in Sect. 6 to the two body problem. Its 
explicit classical solution is presented there. Quantization of the two body problem 
described by nonstandard Schrodinger equation and its solutions, in particular nu- 
merical calculations of energy levels and wave functions in the confinement regime 
are given in Sect. 7. The classical dynamics and quantized solutions in the presence 
of an external constant magnetic field are described in Sect. 8. Sect. 9 reports some 
of our conclusions and describes an outlook for further investigations. 



2 Nonrelativistic Particles and their Covariant Cou- 
pling to Vielbein Fields 



Let us assume that d space coordinates x(t) = (x\(t) . . .Xd{t)) describe a trajectory 
of a ci-dimensional point particle with dynamics invariant under the transformations 
(1.1). The velocities sfcj = ^ then transform under (1.1) as follows: 



2 It has been hypothesized already by Hehl that "torsion should arise in the microphysical realm" 
[11]. Some time ago [12] it was also shown that the gauge theory of gravity with torsion can lead to 
a confining potential. 

3 The preliminary results were presented in [13]. 
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The formula (2.1) can be obtained from the (d+ 1) -dimensional "relativistic" formula 

X,, 

dx' 1 * 



Jfl - (x,x )) 
if we assume that x° = — xq 



x 

—t, i.e 



dx l 



x 



x = 1. 

The preservation of (2.3) in any coordinate frame implies that 



dx' 

dxi 







dx' 
~dt 



Xr\ 



x + a 



(2.2) 



(2.3) 



(2.4) 



in accordance with (1.1). 

It is well-known how to introduce the compensating gauge fields for the trans- 
formations (2.2); namely, we should replace the velocities by the world scalars (i.e. 
scalars under the transformations (2.2)) which are vectors in tangent space 



x" 



(2-5) 



Here the (d + 1) x (d + l)-bein E v ~ transforms as a covariant vector under the local 
transformations (1.1) 



, E _ 8xP n 
K -dx^ Ep 



(2.6) 



and as a global (d+ l)-dimensional vector under the Lorentz rotations in the tangent 
space. The imposition of relations (2.3-4) and their validity in any coordinate frame 
imply that one should choose 



E y - 



M = 



(2.7) 



as then one gets £- = x° = —1. The choice (2.7) we shall call the nonrelativistic 
gauge because it splits the (d + l)-dimensional Lorentz vector into a rf-dimensional 
nonrelativistic vector and a scalar as well as because it implies the Newtonian notion 
of absolute time. One can write (a,b = 1, . . . d) 



( 1, 0, 



\ 



(2.8) 



4 The condition (2.3) can be treated as fixing the general reparametrization x' = x' (x, xo) of 
(d + l)-th coordinate. 
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and introduce the inverse vielbein E u „ as follows: 



( i, o, 

e 1 
V e d 



0\ 



/ 



(2.9) 



where hr i h\ = 6% and e l , e- are independent. 

We replace the "flat" free nonrelativistic Lagrangian by [] 



m 
~2 



£o = —x 2 



m 



C + C mt — — £a 6" ■ 



(2.10) 



Using formula (2.5) for the nonrelativistic covariantized velocity one obtains the 
following Euler-Lagrange equations of motion 



h- p — T- P — T- f ti 
or, using (1.3) and then E % E^q = 0, one gets & = h 3 b £,- — e- 7 and 

ic-h\h\_T%^-h\T? o £. a = 0, 
where the tensors T^, are given by the formulae (see also (1.4)) 



T. a . = h-. - h-. 

L V '%] ' L 3,t 



7^ _ U± _ p a 
1 i0 — '%Q c ,i 



(2.11) 



(2.12) 



(2.13) 



In particular, for d = 1 we have only two zweibein-components Eq = e, Ej = h 
and Zb — 0, — d t h — d x e (see [2]). In d = 2 we have 



hf 



hi, h\ 



hj, h- 2 

and it is easy to see that the term 



(2.14) 



does not vanish. 

One can also write Lagrangian (2.10) for any d as 



fCOV 

^0 



'/??. 



(g ik x l x k + 2g i0 x l + g 00 ) , 



(2.15) 



5 The d-dimensional indices, due to the Euclidean nonrelativistic metric, can be taken equivalently 
as lower or upper indices. 
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where 

g lk = hfhl g i0 = e±hf g 00 = e^e^-l. (2.16) 

We see that one can treat the components of the metric as describing the 
magnetic field coupled to the velocity of the particle^] [18,19], and that the component 
goo describes a potential. The Euler-Lagrange equations of motion now take the form[] 

do 

x l + r jk V x k = -gV J±# + g*i F fl x l + F l [g ij: 9oi , g 00 ] , (2.17) 
where r 4 is the symmetric Levi-Civita connection 

r' i = U^i + p- a -M), (2.18) 

J 2 \ox k oxj oxi J 

the field Fji plays the role of the magnetic field strength (Fji = djA t — diAj, with 
Aj = gj ) and the force F l is given by the formula 

F l = g^E v (2.19) 

where Ej plays the role of electric field strength 

dA dAj 

The most convenient form of the equations of motion depends on the type of field 
action for the vielbein field. The formulation with basic variables £a is the best suited 
for the discussion of the interactions with a gravitational field described by an action 
functional which depends only on vielbeins and on the torsion variables (1.4). 



3 Torsion Lagrangians and Coupled Particle— Torsion 
Field Systems 

Following Moller [20] one can postulate vielbeins as fundamental variables in grav- 
ity theory and treat the metric as a derived quantity. Let us recall here that the 
covariantization of the Dirac equation requires vielbeins. Moreover, tetrads appear 
naturally in the geometric framework as translational gauge fields, providing the so- 
called teleparallel formulation of relativity theory [5-10]. 

Because eq. (2.11-12) involve only tetrads and torsion tensors, it is convenient 
to consider the free gravity actions which depend only on these field variables. In 

6 It is interesting to relate this observation to the Kaluza-Klein description of the electric charge 
in terms of the additional fifth dimension. 

7 The equations (2.16) appeared also in [1], as the nonrelativistic limit of a covariant relativistic 
particle model. 
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order to obtain nonrelativistic actions we can consider the relativistic D = d + 1- 
dimensional torsion actions and set the time component T® u = 0, consistently with 
our nonrelativistic gauge condition (2.7). Because the explicit forms of the torsion 
actions depend strongly on dimensions, we shall consider below, separately, the three 
cases of d = 1, 2 and 3. 

i) d= 1 

It is known that in 1 + 1 dimension the Einstein-Hilbert action is a topological 
invariant and so we can consider the action quadratic in curvature (see e.g. [21]) or 
in torsion [14]. The quadratic torsion action has the form 

sTt =Y\l dtdx ' det E ■ T ^ T ^ { 17 = o, i (3 - 1} 

and reduces to the nonrelativistic field action considered in [2,3] when we observe 
that in the nonrelativistic gauge (2.7) 

detE = h. (3.2) 

Putting T~ v = and writing F = |T~, we get0 

S nrT = — J dtdxh ■ F 2 . (3.3) 

ii) d = 2 

The (2 + l)-dimensional gravity with Hilbert-Einstein action is dynamically trivial 
as outside nonvanishing matter sources the space-time is flat (see e.g. [22]). This fact 
was one of the reasons why the alternative 3-dimensional topological gravity models 
- without torsion [23] and with torsion [16,17] had been considered. In particular, one 
can introduce as a candidate for a D = 2 + 1 gravity action the translational Chern- 
Simons term [16,17] (t = x ; /i, u, p = 0, 1, 2,a,{3 = 0, 1, 2; 7] a p = diag(— 1, 1, 1)) 

S™l = \j d 3 xe^»E%T§- pVaP . (3.4) 
If we pass to the nonrelativistic gauge (2.7) we obtain (a, b = 1, 2) 

StT = {jdtd 2 xe^E±T«-. (3.5) 

The action (3.4) will be used in the next Section to derive and study various properties 
of the planar iV-body interactions. 

As later on we will study in detail the d = 2 case let us describe in some detail 
the planar coupled system of equations describing an interacting particle - torsion 
field system, with torsion fields described by the Abelian Chern-Simons action (3.5). 

8 See formula (18) in [2]. 
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Introducing N trajectories x a (t) = x l a (t) (i = 1, 2; a = 1, . . . N) for N particles and 
the notation 

E^. a (t) = E%x a (t),t), (3.6) 

the free action for N particles in d = 2 dimensions, in the first order formalism, can 
be written as 



S 



(TV) 



N 



(3.7) 



pail ~ m a J dt £, a £, a £ a (hj-a + 6 a / 

a=l 

thus providing the constraint formula 

e = hj. a xi + e% (3.8) 
The action (3.7) can be written in the field form, using the mass density function, as 



SCAT) r M 

part = - y dtd\ y-^e - e- (/^ + e - 



p r (f;fi(t) . . .x N (t)), 



where 



N 



p N (x; xx(t) . ..x N (t)) = rn a 5 i2) (x - x a (t)) . 

a=l 

The field action (3.5) thus takes the form of the field action k = 1, 2) 

S field = \f dtd 2 x (fB* - s jk hjsf) , 

where 



(3.9) 
(3.10) 

(3.11) 

B*- = e jk d s t% £% = d t hl-d k e & . (3.12) 

The fields B- and £ k plays the role of the magnetic and electric fields, respectively, 
with the internal 0(2) index a describing the d = 2 nonrelativistic rotation groupf]. 

We can now derive the coupled equations describing the d = 2 particle-field 
system, described by the action S = S pJt + S field (see (3.7) and (3.9)). The equations 
for the particle trajectories (see (2.11)), having used the notation (3.12), then take 
the form: 

h% e - % B a - e x J - Sf ;a e = 0. (3.13) 

We can now pass to the form (2.11) and use the following explicit form for the 
nonrelativistic dreibein (2.9) 



(3.14) 





( d 





o \ 


E »~d 


6 h 2 6 h 2 1 








\ e-h\ — e-h\ , 


-h\, 


h\ ) 



9 In the vielbein formalism the rotation group (Lorentz in the relativistic case, Euclidean in the 
nonrelativistic gauge) plays the role of an internal symmetry. 
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where d = det (h~j) = h\h- 2 — hr 2 h\. 

The field equations for the dreibein fields then take the form: 

S%(x,t) = -^^nriaSkjia ■ £| • S {2) (x- x a (t)), (3.15a) 
B*(x,t) = ~^2m a &M(x-x a (t)). (3.15b) 

^ a 

We see, therefore, that the dreibein field equations, consistently with the general 
property of the d = 2 + 1 gravity, imply that the space-time is flat outside of the 
matter sources. 

In the next section we shall try to solve the system of equations (3.13) and (3.15a- 
b) by using the technique of singular gauge transformations (see e.g. [24-28]). 
iii) d = 3 

In 3 + 1 dimensions we can write three independent quadratic torsion actions [5- 
10]. It is interesting to observe that there are special linear combinations of quadratic 
torsion terms which provide four-dimensional Hilbert-Einstein action. The detailed 
consideration of coupled nonrelativistic particle - torsion fields systems^] is postponed 
to a future investigation. 



4 Solution of the Field Equations 

4.1 Gauge fixing and residual symmetry 

The set of equations (3.15a-b) can be rewritten as (see (3.5)): 

e ^ Pd ^=*J2^x<>> a 6{x-x a ), (4.1) 
where (3.15a) and (3.15b) are obtained respectively by putting n = 1,2 (// = «) and 

/2 = 0. 

The general solution of (4.1) can be written in the pure gauge form 

E%(x, t)=E , t) + d^(x, t) = d^(x, t) (4.2a) 

where 

E %x, t) = - 53 g $(f - j? o ) (4.2b) 

a 

and 

- A- is an 0(2)-vector valued pair of regular gauge functions, 

10 In our framework, due to the relation (1.4), the local frame fields can be called torsion 
potentials. 
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- <&(#) is a singular gauge function satisfying the following equation (see e.g. [26- 
28]) 

e t3 didj$(x) = 2nS(x ) (4.3) 

thus expressing the singular nature of the first term in (4.2b). As a solution of (4.3) 
we can take 

$(£ ) = arc tan— (4.4) 
x 1 

i.e. 

dk$(x) = SkAhilx |. (4.5) 

The function dk&(x) can be regularized in such a way that it has a well defined limit 
for x — > 0, e.g. we can replace In |x| in (4.5) by [28] 

In \x | — > ln^ (x) := — / rf 2 y In |x — y \e~~ . (4.6) 

en J 

In this case we find 

z) ln (£) (f) ^Mln|f|, 

ii) Yime kl dM e \x) — >0 > 0. (4.7) 

Let us note that the solutions for the fields E^(x, t) with asymptotically nonvan- 
ishing gauge function A- do not solve the Hamilton's variational principle for the field 
action (3.11). The bad asymptotic behaviour for r — > oo of dreibeins E^(x, t) leads to 
the appearence of nonvanishing surface integrals and, in consequence, our E^(x, t) do 
not minimize the action. Such a situation arises also in General Relativity (see e.g. 
[29]). In the following we will argue along similar lines by choosing an appropriate 
asymptotic form for E^ and adding two surface integrals to the action. 

To do this we decompose 

E±(x,t) + E^ (4.8) 

where we require that the new field variables E^(x, t) satisfy 

E%(x,t) -> Oir- 1 ) as r -> oo (4.9) 

and we assume that the asymptotic parts E® s - are given as functions of t only. Then 
from the pure gauge form 

Ep(t) = d^(x,t). (4.10) 

we obtain 

A l (x,t) = x i -a l {t) (4.11) 

giving us 

= -a*(t) =: -v^(t), E™- = Sf, (4.12) 
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where we have chosen the factor in front of x % to be equal to one in order to have, 
asymptotically, the Euclidean metric. 

The choice (4.11) for A* breaks asymptotically the invariance with respect to 
local space translations (1.1). In the formula (4.2a) only the fields E^ transform 
covariantly with respect to those local translations which preserve the asymptotic 
behaviour (4.9). However, as under general coordinate transformations the functions 
A 1 are scalars, the changes 5 A 1 under (1.1) of both the singular and regular parts of 
A 1 must separately vanish and we obtain 

5x i = Sa^t), (4.13) 

where x l and a 1 transform as vectors under rotatons in tangent space. Therefore, as 
a residual symmetry, we obtain translations, local in time, and rigid rotations (see 
also [30-32] for a similar situation in General Relativity). 

Putting (4.8) and (4.11) into the Chern-Simons action (see (3.5) and (3.11)) we 
find (modulo a total time derivative) 

L field = L field — h — h (4-14) 

where 

L feu := Lfieid[E^] (4-15) 
and the integrals are defined by 



I 1= I / d 2 x ■ e ij di Ef(x, t) (4.16) 
I 2 = {fcPxe^diEi (4.17) 

Next we use the Stokes' theorem to rewrite ii j2 as 

h = \ linwoo rv^f^ dip (-Efsin(<p) + £f cos(v?)) , (4.18) 
I 2 = \ lim™ r ; 2 " d V (-El sin(<p) + E% cos(<pj) . (4.19) 

Given the asymptotic behaviour (4.9) the boundary integrals I± t 2 exist but they 
do not vanish, in general, and so only the modified field Lagrangian 

L field = L field + h + h (4.20) 

has well defined functional derivatives with respect to the field E^. 

Due to the asymptotic behaviour (4.9) the variations of E^ and v- are independent 
of each other i.e. the fields E^ and v- appear as new variables. The property that 
in a (2 + l)-dimensional gravity boundary terms give rise to additional degrees of 
freedom has been shown also in [32]. 
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The particle action now takes the form 
s£Z[x,t,E\ = s£Z[x,£,B\ + / dt 



Ee-^-^-Ee 



(4.21a) 



and the modified field action Sgeid — Sfieidt-E 1 ] (s ee (4-20)) is given by S^ R ^ (eq. 3.5) 
but taken as a function of E 

S md [E] = S ( T NR) [E] (4.21b) 

We see that in the action (4.21) there are separated the variables describing the 
"bulk" (Ej^) and asymptotic behaviour (v a ). 

The new variables (v-) describing local gauge degrees of freedom appear as the 
Lagrange multipliers in (4.21a) and are not determined by the EOM (cp. [33]). By 
fixing it as a constant we obtain Galilei invariance as the residual symmetry. 

If we vary E^ we obtain the field equations (4.1) whose solutions should be taken, 
in accordance with (4.9), as 

E%(x,t) = -^0 M ££(t)<lK* -*«(*))■ (4-22) 

By varying S with respect to v- we obtain the constraint 

Ee = 0. (4.23) 

a 

We note that only the presence of the constraint (4.23) yields the correct asymp- 
totic behaviour as well as the correct rotational properties of Eq because 

• from (4.22) we have 

Ei(x,t) -> -A + oir- 1 ), as r^oo (4.24) 

a 

which agrees with (4.9) if (4.23) holds. 

• for a rotation by an angle <p in the 1-2 plane we have 

$ -> $ + p (4.25) 

and so with (4.24) we get 

££-^-*^-E& ( 4 - 26 ) 

which agrees with (2.6) only if (4.23) holds. 
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4.2 Conservation Laws 



Due to the constraint (4.23) the total momentum of the system vanish. To see this 
let us note that after gauge fixing our Lagrangian L is invariant with respect to space 
translations local in time, and rigid rotations. According to Noether's theorem, the 
invariance of L with respect to 5x leads to a conserved quantity C[5x\ (see [34], eq. 
(6.29c) for the corresponding case of the Chern-Simons electrodynamics): 

CM = E^ • 5 < + \ j d 2 xB^Ef5x k = C part + C field . (4.27) 



account the definition of the canonical particle momentum = j^w we get 



Inserting into (4.27) the expression (3.15b) for B- and (4.22) for E^ - and taking into 

dL 

n . a 

cm = Ee^ (4.28) 



Taking Sx as describing respectively space translations (&c| = 5a-) and rotations 
(&c| = e — x| • 5a) and denoting the corresponding conserved quantities by P i and 
J, we obtain due to (3.15b) and (4.22) for P { the formula 

P l =^p art = EC (4-29) 

a 

and for J we get 

J y ' 6jj Xj jQ , C,j : a- (4.30) 

a 

Observe that J is different from the total canonical particle angular momentum J par t- 
The field contribution J field-, which is obtained from (3.15b) (4.22) and (4.23), namely, 

j field = -^E£e« (4.3i) 

is conserved separately because it is proportional to the ^-particle Hamiltonian to 
be given in section 5.2. 

The fact that J ^ J pa rt will play an important role in the quantization of the 
system (see Sect. 7). 

In section 3 we showed that the residual symmetry contain the translations local 
in time. The conserved quantities P l , which vanish according to (4.23) and (4.29) in 
the physical phase space are the generators of this symmetry in an extended phase 
space. The same result has been obtained for (2 + 1) -dimensional gravity [35] and by 
one of the present authors (PSC) in the one-dimensional case [2], [3]. For the similar 
case of local (in time) rotations it has been shown that the corresponding rotation 
generator is a first class constraint which vanish in the physical part of the phase 
space (cp. [36], [37]). 
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5 Classical Particle Dynamics (N Body Problem) 
5.1 General Properties of the Solutions 

Let us look at the classical equations of motion and their solutions. The equations to 
solve are then 

and 

e • E% + e • F% a x^ a = (5.2) 

for the particle coordinates {x a (t)} after the insertion of the solution (4.2) for E% 
into these equations. 

However, from (4.1) we obtain for F^ v 

F% = -^wEef^^-^)- (5.3) 

Thus the second term in (5.2), due to the antisymmetry of the e tensor, can be 
rewritten as 

e • F%^ a = ~^i P E $ ^x^5(x a - x p ) (5.4) 

which is infinite for coinciding particle positions and vanishes otherwise. Therefore 
our configuration space contains only noncoinciding particle positions. 
From eq. (5.2) and (5.4) we have 

e • E% = 0, (5.5) 

which leads, for points in the configuration space where the metric is non-degenerate, 
to 

e = 0. (5.6) 
Using this result we obtain for given by (1.3) from (4.2) and (4.10-12) 

e = ^-« a -zES • ma*/,)) , (5.7) 



X a0 '■— X a X 0- (5-8) 



where we have defined 

For the particular case of a particle motion on a line we find, from the last two 
expressions, that 

x% = (5.9) 

i.e. the free motion. 
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This is also easily seen from (see (1.3c) in [26]) 



xl p ■ d a $(x a/3 ) = x aP A = (5.10) 

as x a/ 3 and XqlQ £1X6 parallel for a motion along a line. 

We should say also a few words about the degenerate case i.e. when (5.6) is not 
necessarily satisfied. 

Then we have for from (4.2) and (4.10-12) 



e = ^ - v a - - - e e P ■ - t-e - $(o)) , (5.H) 

where $(0) has to be defined by an appropriate regularization procedure. But for any 
such procedure the last term in the last formula breaks rotational invariance. Thus, 
although £f (x, t) is rotationally invariant this is not the case for 

E% = Bm (5.12) 

if (5.6) does not hold. 

Hence the degenerate case can be considered as unphysical and further will not 
be discussed. 



5.2 Reduced Hamiltonian/Lagrangian for the Non— Degenerate 
Case 

Applying the Legendre transformation to the total Lagrangian given by (4.21) we 
find the following Hamiltonian 

H = -J dx 2 E%- feidiEf + 

+|Eee + ^-Ee (5.i3) 

a a 

The gauge field Eq is a Lagrange multiplier field, whose variation gives the con- 
straint 

e^d t Ej(x,t) = (5.14) 

which corresponds to the fx = term in (4.1). 

When if we put this expression into the Hamiltonian we find 

H {N) = ^Eee + ^-Ee (5-15) 

" a a 
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This is a free Hamiltonian with a constraint and the dynamics is contained in the 
non-trivial symplectic structure. 

In order to complete the Legendre transformation we have to express the auxiliary 
variable in terms of the canonical variables {x^, p|}. This can be done in a unique 
way iff 

' dp a 



det [ ^ | ^ 0, (5.16) 



where 



Pa ' Br.i ' ai 



dx l a 

The condition (5.6) defining the non-degenerate metric now arises as a consequence 
of the EOM. When we apply an inverse Legendre transformation to (5.15) we get 

Lred = X>|-*i - H red . (5.18) 

a 

Varying S re d with respect to x l a we obtain using (5.17) 

which leads only to the trivial solution 

t = 0, (5.20) 

if (5.16) holds. We conclude therefore that the existence of the Legendre transforma- 
tion implies a non-degenerate metric. 

Varying S re( i with respect to v- and £g we obtain, respectively, the constraints 

Ee = 0, (5.21) 

a 

in accordance with the relation (4.23) obtained in geometric way by suitable gauge 
fixing and 

e = i£ - - ^ E (5 • • «*)) • (5.22) 

Alternatively, we may consider the Hamiltonian (5.15) in the {x^, £§}-space. In such 
a case the free Hamiltonian (5.15) will be endowed with the nontrivial symplectic 
structure. 



11 Compare with [34] where a similar discussion of the Chern-Simons electrodynamics for point 
particles is given. 
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5.3 Transformation to flat coordinates 

Let us recall from Sect. 2 that the space parts of our dreibeins Ef determine the 
metric in the two-dimensional space 

gij (x,t) := (Ef(x,t) ■ Ef(x,t)) (5.23) 

with the line element given by 

ds 2 = gij dx i dx j . (5.24) 

But according to (4.2) the Ef- are singular due to the formula (we put A- = y-) 

Ef = d lV ^x,t) (5.25) 

with the singularities located at the particle positions {x^(t)}, (a = 1, . . . N). We see 
that our metric is locally flatQ and we have 

ds 2 = {dy^ ■ dy^) (5.26) 

in R 2 -{x±(t)},(a = l,...N). 

We can now show that the variables {y%(i)} are canonically conjugate to the 
covariantized velocities To show this we consider F({£%,x%}) defined as 

follows: 

f := E £ • < - h E (ef • e) (5.27) 

F is a generating function for the canonical transformation 

{*5,Pg} - (5-28) 
at the points in phase space for which 

„et (H-) + 0. (5.29) 
Indeed, the relation (5.17) takes the form 

dF 

Pa,i = TT— (5-30) 

and from (5.25), (4.2) and (4.8-11) it follows 

dF 

y a ,i = (5-31) 
Thus we conclude that the motion is free in the part of the phase space {?/§, £§}^£| =0 

a 

restricted by the condition (5.29) and with substracted coinciding particle positions. 
12 Compare with (2+1) dimensional gravity (cp. [22]) 
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6 Two Body Problem 

6.1 Symplectic Structure 

We define: 

e:=~(£f-£f), x^.= xf-xt := i(p? - pf) . (6.1) 

Then using the constraint (5.21) we get from (5.15) 

H m = ^ . ( 6 .2) 

and 

Pi = 6 + ^(e-e-)^$(f). (6.3) 

The Hamiltonian equations take the form 

ii= ™ = J?.m, (6 . 4a ) 



cfy? 2 \ dp 1 



A = -£=-*(r-£V (6.4b) 



Using (|6.3| ) we have 



? -f? = — jt 1 • < 6 - 5a > 

<* i + A (e . s<i $) 

e .|T = a MW* . (6 . 5b) 
8x1 4n i + ^(C ■ a.*) 



Taking the time derivative of (|6.3|) and using (6.4-5) we get 



£i + ^d&(x){ j £ j = 0. (6.6) 



Let us now illustrate for N = 2 the procedure leading to the free motion (5.20). 
Instead of the canonical variables (xi,pi) we can use the variables (a;*, Then the 
Lagrangian obtained from the Hamiltonian (6.2) would have had the form: 



Pi(£i,Xi) ■ ±i- H 
A 

All 



+ ) i, _ . (6.7) 



18 



The variation with respect to & is given by the expression 

6 = , (6-8) 

1 - - 

which is equivalent to the Hamiltonian eq. (6.4a) with the insertion of (6.5a). Insert- 
ing (6.8) in (6.7) we get 

L=\—? ■ (6.9) 

1 - s (i-m 

In particular if we observe that 

we see that if the velocities are expressible in terms of the canonical variables then 
from (6.6) it follows that 

6 = 0. (6.11) 

Using the Hamiltonian (6.2) we get for the pair of noncanonical variables the Hamilton 
equations (6.4a) in the form xi = {x t ,H} as well as (6.11) given by & = {£i,H}, 
provided that we assume the following nonstandard symplectic structure: 

{•'•'••'•'} U-U 0. (6.12) 

{•<•'•;,} 2 4^ • (6-13) 

1 + 2^ ( W&) 

It is easy to check that the Poisson brackets (6.12-6.13) satisfy the Jacobi identity. 
6.2 Two Conserved Angular Momenta 

For our system we have two conserved scalar angular momenta (cp. section 4.2): 

i) If we define (in d — 2 a A b = e^d'ty is a scalar) 

l:=xAp=xA£ + -^-H (6.14) 

47T 

we find that / is conserved 

ii) Second conserved angular momenta I is the following 
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Z:=£A£ (6.15) 

because 

i = £ a|* = 0. (6.16) 

Using the relation 

e-<9 a $ = -f AVlnr = ^, (6.17) 
where r := \x\, we see that (6.8) can be rewritten as: 

2^ 

x- = — V- (6.18) 

1 ~T 27rr 2 

— # — 

Note that if x(0) is parallel to £ we have / = and, in consequence, a free motion 
on a line. 

For AZ < it is convenient to introduce the quantity 

rl := § (6 . 19) 

We see from (6.18) that the relative two-body problem separates in this case into a 
motion within the interior region given by 

r < r (6.20a) 

and a motion in the exterior region given by 

r > r . (6.20b) 
6.3 Structure of the classical phase-space M 

Let us denote by M the phase space for the canonical variables (x,p). First of all, let 
us observe that 

M ^ R 2 <g> R 2 . 

To show this we start with (6.3) rewritten as 

f=p-AfTv0. (6.21) 

Squaring it we get 

( \ H (2nr l\\ 2 „ /2yrr A 2 
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where we have used the definition (6.14) of I. Thus we must have 

2 



o 2%r I 



< 



(6.23) 



and so (with I = rp sin ip) 



r > (1 — sin ip), for A > 0, 

p(l + sm<^) for A < 0, 



r > — 

' — 2tt 



(6.24) 



where the equality sign holds for r = r . Therefore the map (x, £) — > (x, p) with 
(x, G -R 2 (8) i? 2 leads to (f,p)GM^fl 2 ® # 2 . 

From (6.22) we may look for the Hamiltonian H as a function of the canonical 
variables (x, p) 

4nr 



H 



A 



2vrr I 
A r 



/27rr P 



p 



(6.25) 



with the — (+) sign in front of the second term in (6.25) for (^y- + > (< 0). Let 
us note, however, that due to formula (6.14), the RHS of (6.25) depends through / 
also on the energy H. Writing (6.22) as 



we get 



H 



'I 2nr 



47rr 



I A r 47rr / 



2vrr A 

— + ; =■= 



\ 



p 2 + 



2yrr / 



(6.26) 



(6.27) 



The formulae (6.25) and (6.27) can be used alternatively in the quantization proce- 
dure, depending on the choice of the eigenstates of angular momenta (/ in the case of 
eq. (6.25) and 7 in the case of eq. (6.27)). 



6.4 Classical Trajectory 

To find the classical trajectories we can make the ansatz 

x(t) = a(t)i+ /3(t)£ ± , 



(6.28) 



where we have written x for x- and similarly to £ and defined £± as a unit vector 
perpendicular to £. Then 

£ A £l = \£\ = £. (6.29) 
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As £ is conserved, the conservation of / tells us that j3 is constant and so ( |6.28p 
becomes 

x(t) = a(t)g + /3g± (6.30) 

with a constant (3. So the problem has been reduced to finding a(t). To do this we 
choose a frame in which £ = £ e x . Then £j_ = e y and we see that 

Six) = arctan f — r-r- ) . (6.31) 

Next, as in the previous section, we introduce y 

y = x + ^i<P(x), (6-32) 

which is canonically conjugate to £, and find that using the last two formulae we see 
that a(t) satisfies 

2 ^ + y(0) = a(t)i+ (3g ± + arctan (^7) , (6-33) 

ie. we arive at a fixed point equation for a(t) 

a(t) = — — arctan [ — ^r- ) + 2t + c, (6.34) 
2tt \a(t)£J 

with c determined by the initial conditions. 

To discuss the conditions on the existence of a(t) which solve ( |6.34j) we rescale it 
by introducing 

Then after a time translation t + c — > £ our equation ( |6.34j ) becomes 

9 = f(g) + r, (6.36) 

where 

f(g) = k arctan ^-^ . (6.37) 

To discuss the solvability of (6.36) it is convenient to redefine / to be given by 

f(g) = —k arctan g (6.38) 

and perform a further time translation so that the equation for g is still of the form 
(6.36). Then putting 

g = v + r (6.39) 
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we find that (6.36) becomes 

[■v+t 

= —k I dx 

1 + x 



y+T 1 

k I dx =: g T (u). (6.40) 



Note that this implies that 

v( — (Yll = k 

2 



//(-oo) = k^- (6.41) 



and 

u(oo) = ~ k \ ( 6 - 42 ) 

independently of the value of k. 
Note further that 

7T» = - fc T-^ si ( 6 - 43 ) 

dv y ' 1 + (u + r) 2 v ' 

Let us consider first the case of k > — 1. 
We rewrite (6.40) as 

G T {y) := v - g T = 0. (6.44) 

Then from (6.40) we have 

G r (±oo) = ±oo. (6.45) 
Differentiating (6.26) with respect to r we find 

^ = - mr^TTTi ( 6 - 46 ) 

(i/(r) + r) z + 1 + 

From (6.43) we conclude that G T {y) is a monotonically increasing function of v. 
Given (6.45) and (6.46) this implies that there is a unique and smooth solution v{j) 
of (6.44) for each t e R 1 . 

For k < — 1 we have a problem as v develops a discontinuity. This can be seen 
by plotting the various terms in (|6.36|) . Thus we see that for k > — 1 we have a well 
defined g{r) while for k < — 1 the function p(r) "jumps" at some value of r showing 
that the x space is not the physical configuration space in this case. 

Let us discuss this in more detail. 

Eq. (6.46) shows the smoothness of v(t) for all values of k > —1. For k < — 1, 
however, u(t) has an infinite slope at r = ±r , where 



To :-- 



\k - 1 ± A; / dx -. (6.47) 

Jo 1 + x 2 



But due to the symmetry u(—t) = —v{t) it is sufficient to consider the case of 
t = tq. To do this we note that if we start at r = — oo we have a smooth function 
v{t) for r < tq. Close to r = r we have 

lim z>(r - e) = +oo (6.48) 
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with 

i/(ro - e) = -r - y/\k\ - 1 + e (6.49) 
thus showing that u{t) jumps at r = r to 

z/(r + e) > i/(r -e). (6.50) 

However, the states z/(r ± e) are physically equivalent; i.e. they have to be 
identified in the configuration space C p h ys - (6.33) shows that C p h ys is determined by 
l/-space defined by (6.32). To see this we recall from section 5.3 that the map x — > y 
is part of a canonical transformation iff 

(d 2 F \ 
a**) * (6 51) 

which for iV = 2 is equivalent to 

l + ^C^^O (6.52) 

and so 

{v + r) 2 + 1 + k ^ 0. (6.53) 

But for the whole trajectory i/(r) we have (with the identification v(tq — e) ~ 
K^b + e)) 

[v + r) 2 + 1 + k > 0, (6.54) 
which, due to (6.17) is equivalent to 

r 2 > r 2 for XI < 0. (6.55) 

Thus we have a one-to-one correspondence between the y-sp&ce and C p h ys for 
r > r . 

Within the interior of the region r < r$ we obtain from (6.40) a solution valid 
only for a finite time interval. Such a situation is well known in General Relativity 
(cp. [38]) and the literature cited there). 

From the physical point of view the case k < — 1 is the most interesting as then 
the configuration space consists of two non-communicationg regions; the exterior 
(r > r ) and the interior (r < r ). Thus, in the next section we will concentrate on 
the quantization of this case. 

7 Quantization of the Two— Body Problem 

7.1 Nonstandard Schrodinger Equation 

The relation (6.3), after the substitution of (6.2), takes the form (H = H^>): 
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A 



f = P 

Squaring it and using again (6.2) we obtain 

I 2 



H = p 2 



+ 



(7.1) 



(7.2a) 



where we have used the definition (6.14) of 
Let us note that 

1 = 1-— H 
Air 

i.e. (7.2a) gives us a quadratic equation for H (cp. section 6.3). 

We quantize the problem by considering a Schrodinger-like equation 



ih — — — = Hifj(x,t) 



I 2 

P 2 -- + -l [ 



(7.2b) 



(7.3) 



in which the operators H and p are defined by the usual quantization rules 



Pi = -di. 

i 



(7.4) 



We see that the equation (7.3) describes a nonstandard form of a time dependent 
Schrodinger equation, with its right hand side containing both the first and second 
time derivatives (entering through 7) . 

For the stationary case, i.e. when ^(x,t) = ^/e{x)^ k ' we can use the angular- 
momentum basis and put 

y E , m = f E , m {r) ^ (7-5) 

where m is an integer, and find that fE, m satisfies a nonstandard time independent 
Schrodinger equation 



-h 2 ( d 2 + -d, 



— i s 

m 



- E 



fE, m (r) = 0, 



(7.6) 



where in consistency with (7.2b) we have defined 



km 



km — —E 

47T 



(7.7) 



i. e. hrh is an eigenvalue of 7. 

A characteristic feature of the Schrodinger equation (7.6) is the appearance of the 
noncanonical angular momentum Tim with m not being an integer, i.e. our two- 
body system carries a fractional orbital angular momentum (see the discussion in [27] 
or [39]). It has been shown in section 4.2 that J = I is equal to the total angular 
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momentum of the particle + field system with a nonvanishing angular momentum 
of the fields. We note that a form of (7.7) shows a great similarity bteween our 
two-particle state and the gravitational anyon of Cho et al. [40] as an energy-spin 
composite. 

In the following we discuss only the most interesting case of XI < 0. 

Now the appropriate boundary conditions correspond to the requirement that 
fE,m{ r ) is nonzero in either the interior region (r < r ) or in the exterior region 
(r > r ). Thus our boundary condition is 

f E ,m(ro) = (7.8) 
The general solution of (7.6) is given by 




Mr) = 4 V r > ( 7 - 9 ) 



where Z w is an appropriate Bessel function of order m (or a superposition of such 
functions). 



7.2 Interior Solutions (r < ro) 

The only Bessel functions not blowing up as r — > 0, i.e. giving a finite probability 
f 2 r dr as r — > 0, and possessing positive zeros r are those of the first kind of order 
|m| with E > 0. 

Then the possible eigenvalues E n (m) are determined by 



J 7) 



E ( Ti\Xm\ 



2tt 



(7.10) 



with m given by ( |7. 7\ ) . 

Let us look in more detail at the case rh > 0, A < 0. To simplify ( [7. 10| ) we define 



\X\E 
2-nh ' 



(7.11) 



Thus (|7.10p takes the form 



0. 



(7.12) 



As Jm, for fixed rh > 0, has an infinite number of positive zeros, which we denote 
by y n (m), n = 1,2.. we see that due to (|7.7|), the eigenvalues e n (m) we are looking 
for are the positive fixed points of the equation 



f n (m+ -e), 



(7.13) 
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where we have defined 



f n (m) = -y 2 n (m). (7.14) 
m 



The existence of positive fixed points e of ( [7.1 3D may be shown by using appropriate 



bounds for f n or y n . For example for n = 1 and m > we have [41] 

4 

m(m + 2) < mfi(m) < ~(m + l)(m + 5) (7-15) 

3 

and hence we have, with m = m + and for m e iV, 



2(m + 2) < ei(m) < m + 12 + (m + L2) 2 + 8(m + l)(m + 5) 2 , (7.16) 

i.e. we have proved the existence of €\(m) > for m = 1,2.. and give crude bounds 
for it. A better estimate for large m can be obtained by using the asymptotic formula 
[41] 

yi(m) = m + 1.855757m 3 + 0(m a) 



valid for large m leading to 

e x {m) = 2m + 9.35243m^ + 0(rn~%) (7.17) 
valid for large positive m. 



To get more insight the equation ( 7.13 ) has to be solved numerically. To do this 
we have to determine the behaviour of zeros of the Bessel function as a function 
of k. Then having determined this dependence we can find the values of energy by 
the secant method. 

Luckily there are many computer programs to determine zeros of Bessel functions 
and in our work we have used the Maple program to perform this task. As the 
dependence of each zero is almost linear the numerical procedure of solving ( |7.13| ) 
converges rapidly. 

In fig.l we present the values of energy for e < 300 as a function of m. The 
plot looks like several curves; the lowest values correspond to first zeros (ie n — 1), 
the next ones to second zeros ie n = 2 etc. The points lie so close that the figure 
may appear as a set of lines while, in reality, we have here sets of points. The points 
appear to be (almost) equally spaced on each "curve" - this is due to the approximate 
linearity of the positions of zeros of Bessel functions as a function of m. To check 
our values of energies we have also solved ( |7.13|) differently; we approximated the 
positions of zeros of the Bessel functions by a linear function and solved the resultant 
equations for e. The obtained results were very similar to what is shown in fig. 1 
thus giving us confidence in our results. 

Our results show that, for each value of m, there is a whole tower of values of 
e corresponding to different zeros of the Bessel functions. The values of e increase, 
approximately linearly, as we take higher zeros (ie y n for larger n). The dependence 
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on m is only slightly more complicated; for each order of the zero there is a value of 
m for which the energy is minimal and as we move away from this value the energy 
grows, approximately, linearly. As n increases the minimal values of m increase, 
again, approximately linearly. 




Figure 1: Energy as a function of m 

Note that our numerical values of e are consistent with the asymptotic results 
mentioned before. Note also that for m < and A > the corresponding energy 
levels are obtained by changing the sign of m. 

Summarising, we see that in the interior region r < r , where classical solutions are 
only possible for a finite time interval, we find quantum solutions which correspond 
to discrete bound states determined by the boundary condition at r = r^. Thus we 
see that this boundary condition defines a geometric "bag" for the quantum state. 
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7.3 Exterior Solutions (r > 7*0) 

First of all, let us note that there are no bound states solutions of (7.6) for r > r as 
the only square integrable Bessel functions in [r , 00) are the modified Bessel functions 
of the third kind, which have neither positive nor pure imaginary zeros (see e.g. [41]). 

Scattering solutions are given by a superposition of Bessel functions of the first 
and second kind 



f E , m (r) = A m (E) J<n \^rj + B m {E) Y m \^rj (7.18) 

with the ratio ^ determined by the boundary condition (7.8). These solutions 
describe scattering on an obstruction of radius r , which is dynamically determined. 

8 Two Charged Gravitationally Interacting Parti- 
cles in a Magnetic Field 

In this section we consider the motion considered in Sect. 6 and 7 of two particles, of 
equal electric charge e, under the influence of an additional static uniform magnetic 
field B perpendicular to the plane of motion. We assume, for simplicity, that eB > 0. 

8.1 Classical Dynamics 

To have the Lagrangian L, describing the relative motion of two charged particles 
in an additional constant magnetic field B we add to the Lagrangian considered in 
Section 6 the term ^ 

Lmagn ^ 6 B €ij Xi Xj . ( 8 . 1 ) 

We obtain 

L = (& + — fd^)^ - e 2 + -eBt l] x l x j (8.2) 
and therefore, in accordance with (6.2) we have 

H = i 2 . (8.3) 
The equations of motion for x iy derived from (8.2), are the same as before, ie 



where the covariant velocities & for B ^ are not conserved. However, as & satisfy 

(8.5) 



x 
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we see (from (8.4-5)) that the conserved quantities are now 

eB 

& := & Y €ijXj ' ( 8 ' 6 ' ) 

From the equations of motion (8.4-5) and the Hamiltonian (8.3) we derive the 
following, nonstandard, symplectic structure for the noncanonical set of variables 

(x,q) or (x,i~), respectively: 



{xi, xj} = (8.7a) 

(8.7b) 



1 + n=?iOl 







« = TT^ (8 ' 7c) 

but 

=-^eSey. (8.7d) 

It is easy to check that the Poisson brackets (8.7a-d) satisfy the Jacobi identities. 
With the canonical momentum pi derived from (8.2) 

A cB 
Pi = & + H di4> -^tijXj (8.8) 

we find that the relation between the two conserved angular momenta 1,1 is not 
modified 

l: = xAp = I + —H; (8.9) 

47T 

however, / is given by 

l:= x A|* + |eBr 2 . (8.10) 

We note that due to (8.9-10) the form of the denominator in the equation of motion 
(8.4-5) is different when compared with the B — case: 

1 + s M *=( 1 -s efl ) + 2^- (8 ' n) 

Let us consider just the nonconfined motion, ie (8.11) should be nonvanishing for 
all < r < oo. When B = this would correspond to the case of XI > - ie "the 
less interesting case" of the previous section. 

Now, however, we have to distinguish two subcases: 

Either 

A : 1 - — eB > 0, (8.12a) 
8ir 
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which allows both signs of A, 

and XI > (8.12b) 

or 

B : 1 - — eB < 0, (8.13a) 

8tt v 1 

which allows only A > 

and I < 0. (8.13b) 

Note that when B ^ we can have a nonconfined motion also for XI < 0. 
In analogy to the free field case we introduce for 

A/ < (8.14) 



the boundary r between the interior (r < r ) and the exterior (r > r ) region defined 
by 



2 1 



XI 



(8.15) 



When we consider the confined motion we consider the motion in the interior region 
(r < r ). In this case we have to consider two further subclasses: 
Either 

C : 1-^eB > 0, (8.16a) 



which allows both signs of A, 



or 



8tt 

and XI < (8.16b) 

D : 1 - — eB < 0, (8.17a) 
8n 

which allows only A > 

and I > 0. (8.17b) 

Note that when 5^0we can have a confined motion also for XI > 0. 

Comparing all four subcases we conclude that we obtain confinement <-> noncon- 
finement transitions by tuning the strength of the magnetic field B. To be more 
specific let us consider eg X > and / > 0. Then we obtain, for eB increasing from 
low values towards a transition from case A to case D ie a nonconfinement — > 
confinement transition. This transition is continuous because we have r = oo at the 
transition point. 

In order to quantize our system we have to generalize for B^0 the relation (7.2a). 
We start with (8.8) rewritten as 

& = Pi - ^-Hdi4> + ^-etjXj, (8.18) 
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square it and then using (8.3) and (8.9-10) we find for the Hamiltonian H 

H = p 2 - 4 + 4 + 4( eS ) V " \ eBL ( 8 - 19 ) 
r z r l 16 2 

8.2 Quantization 

Following the method presented in Section 7.1 we obtain from (8.19) the following 
nonstandard Schrodinger equation for the radial wave function f E<m 

f E ,m(r) = 0, (8.20) 

which generalizes (7.6). Like in the previous case, the noninteger eigenvalue fh is 
related to the integer m and energy E by the eq. (7.7). 



U H- h - £ U ) r > - E - §eZM 



8.2.1 Nonconfined motion 

The eigenvalue problem (8.20) is identical to the two anyon problem in a constant B 
field with the statistics parameter being proportional to the energy value E. So, for 
the energy levels E n)Tn we obtain ([42-44]) 

1 X Ti 

E n , m = %eB(n + -\m\ + -) - -eBm, (8.21) 

with n = 0, 1,2, .. Thus, according to the two cases A and B defined above and for 
different signs of A we have to distinguish between three cases: 

Case A with A > 

Due to (8.12b) and (7.7) we have 

fh = m - -\ E n . m > 0. (8.22) 



E n<m = heB(n+^) (8.23) 



Therefore we get from (8.21) 



corresponding to the case II in i43fl . 
Combining (8.22) and (8.23) gives us 

XeB ( 1\ . 
m > -to [ n + 2) ' (8 ' 24) 

where, due to (8.12a), the product of the coupling strength A and of the field strength 
B is bounded 

XeB < 8n. (8.25) 
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Case B with A > 

Due to (8.13b) and (7.7) we have 



Thus, from (8.21) we obtain 



m = m ^- E nim < 0. (8.26) 



E„.n, Th-D ( n + \fh\ + - ) . (8.27) 



which corresponds to the case I in [43 



Inserting (8.26) into (8.27) we obtain 

heB 



E n,m = Y~EEE [ n ~ m + 2 J ■ (8 ' 28) 

Note that in order to satisfy (8.13a) the product of the coupling strength A and 
of the field strength B must be above the mininum value 

XeB > 8tt. (8.29) 

Therefore we obtain from (8.26) and (8.28) 

m> ^( n + \\ (8.30) 



4tt V 2 

Thus we see that E n>rn is bounded from below by 

E n>m > heB(n + -). (8.31) 

Case A with A < 

The condition (8.12a) is now fulfilled automatically. From (8.12b) we obtain (8.26) 
and therefore we have the result identical to (8.28). 
Now we obtain 

m<-^(n+i) (8.32) 
and E n ^ m is bounded from below again by 

E n , m > TieB{n+ 1 -). (8.33) 
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8.2.2 Confined motion 



Now we have to consider a nonstandard Schrodinger equation (8.20) in the interior 
region (r < r ) with the boundary condition as in the previous case 

f E , m (r ) = 0. (8.34) 

The regular solution of (8.20) is given, up to a normalisation factor, by [45] 



E,m(r 



m + 1 



7, 1 + \m\; (3r z 



(8.35) 



with 



and 



1 Ah 



E 



heB 



(8.36a) 
(8.36b) 



and where (p(a, b; z) denotes the confluent hypergeometric function. 

From (8.15), (8.34) and (8.35) we conclude that our energy levels E n ^ m are deter- 
mined by the roots of the equation, which only be solved numerically, 



m + 1 



E „ ._. eB 
,1+ m; — 

8"7T 



V 2 heB 

where m is a function of E as given by (7.7). 
The zero B field limit 

From the well known relation [45] 



Am 



0, 



(8.37) 



lim <P{a,b;--) = T{b)x^ {1 - b) J b ^{2^) 



(8.38) 



we obtain, from (8.35) and (8.37), in the vanishing B limit, the considered in section 
7, respectively, the wave function and the eigenvalue condition. 
The high B field limit 

Without going into the numerics we can conclude from (8.37) that for eB — > oo 
the energy levels increase with the increase of eB, at least, linearly. We prove this 
statement by assuming the contrary. Then one has at the l.h.s. of (8.37), with m > 0, 
due to (8.16-17), <j>{\, 1 + m; m) for which we have the inequality 



-,l + m;m) > 1 for m > 



5.39) 



in contradistinction to (8.37). 
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Confinement — > Nonconfinement transitions 

According to (8.15) confinement — > nonconfinement transitions occur if r — > oo, 
ze if 

1 - — eB -> 0. (8.40) 
8tt V ; 

But from the asymptotic behaviour of the confluent hypergeometric function [45] 

<f>(a,b;x) -^oo ^y^" 6 (8.41) 

we conclude that the r.h.s. vanishes only if a = 0,-1,-2,..., ie for (8.37), if the 
energy levels are given by (8.21) as required. 



9 Conclusions and Outlook 

The aim of this paper has been to consider, for two-dimensional nonrelativistic par- 
ticles, a new interaction scheme, generated by the coupling to the topological torsion 
Lagrangian. We have shown that only in D = 2 + 1 one can write the torsion La- 
grangian in the form of a bilinear Chern-Simons term, described by the translational 
gauge fields (threebeins) multiplied by their Abelian field strength (torsion fields). 
By a suitable choice of the reparametrization gauge, we have then determined the 
solutions of the two-body dynamics with fractional angular momentum eigenvalues 
corresponding to trajectories confined to finite regions of d = 2 space. The quanti- 
zation problem was described by a new type of Schrodinger equation, with a second 
order time derivative, which in the stationary case gave us a nonlinear energy eigen- 
value problem, describing infinite sequence of bound states. 

The eigenvalues of energy were determined numerically and we showed that the 
spectrum was descrete and characterised by two integers; one of them corresponding to 
the rotational integer quantum number m and the other, another integer - n, described 
which zero of the Bessel function the state corresponded to. The dependence on both 
of these integer-value parameters was approximately linear - as we discussed this in 
section 7. 

We have also noted that due to the topological form of the free gravitational field 
action the case of two dimensions is exceptional. In d = 3 one can also look for a 
modification of the standard gravitational interactions (e.g. of the Newton potential 
in the lowest order static approximation) by adding to the Hilbert-Einstein action 
appropriate bilinear torsion terms (see e.g. [7]). However, the confinement due to 
the presence of torsion in D = 3 + 1 gravity leads in the three-dimensional space to 
an additional potential term of a harmonic form (see e.g. [8]). Such a mechanism 
of confinement is different from our proposal, where the two-particle Lagrangian be- 
comes singular at the boundary r of a dynamically determined compact region of 
space. Equivalently, confinement is obtained from the singularity of the nonstandard 
symplectic structure within the Hamiltonian formalism. 



35 



Our considerations have thus shown us [2,3,13] that the geometric bag formation^] 
based on our dynamical assumptions are possible in d = 1 and d = 2 spaces, but 
it is not clear how to obtain analogous solutions, to be generated by nonstandard 
gravitational interactions, in the d = 3 case, i.e. in the standard physical space- 
time. It is worth pointing out, however, that our geometric bag solutions carry some 
resemblance to the fields describing black holes (see e.g. the dynamical division of 
space into two disconnected domains - interior and exterior). At present we can 
only hope that all such consequences for d > 2 will be further clarified in the future 
development of our approach. 
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